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Abstract 

The properties of string networks at scales well below the horizon are poorly under- 
stood, but they enter critically into many observables. We argue that in some regimes, 
stretching will be the only relevant process governing the evolution. In this case, the 
string two-point function is determined up to normalization: the fractal dimension 
approaches one at short distance, but the rate of approach is characterized by an ex- 
ponent that plays an essential role in network properties. The smoothness at short 
distance implies, for example, that cosmic string lensing images are little distorted. 
We then add in loop production as a perturbation and find that it diverges at small 
scales. This need not invalidate the stretching model, since the loop production occurs 
in localized regions, but it implies a complicated fragmentation process. Our ability 
to model this process is limited, but we argue that loop production peaks a few orders 
of magnitude below the horizon scale, without the inclusion of gravitational radiation. 
We find agreement with some features of simulations, and interesting discrepancies 
that must be resolved by future work. 



1 Introduction 



The evolution of cosmic string networks is a challenging problem. The need to consider large 
ratios of length and time scales makes a complete numerical analysis impossible, while the 
nonlinearity of the process defeats a purely analytic treatment. Thus a full understanding 
of most of the signatures of cosmic strings will require a careful combination of analytic and 
numerical approaches. This is needed both to establish precisely the current bounds on the 
dimensionless cosmic string tension Gfi, and to anticipate what will be the most sensitive 
future measurements. Also, if cosmic strings are one day discovered, a precise understanding 
of the network properties will be needed in order to distinguish different microscopic models. 

On scales close to the horizon size, the networks are reasonably well understood from 
simulations DEiQ In particular, there are a few dozen long strings crossing any horizon 
volume. On shorter scales, however, the situation is far less clear. There is a nontrivial 
short distance structure on the strings, which arises because the intercommutation process 
produces kinks [6]. There have been many previous analytic and numerical studies of this, 
but there is no general consensus as to its nature. 

In particular, the size at which typical loops form is uncertain to tens of orders of mag- 
nitude. One widely-held assumption has been that gravitational radiation is necessary in 
order for cosmic string networks to scale [7j, and that it determines the size of loops. If so, 
the loops will be parametrically smaller than the horizon scale, by a power of Gn [Sj (even 
shorter scales have been suggested [9]). More recent work appears to be converging on loop 
sizes at a fixed fraction of the horizon scale, but even here there are significant differences. 
Ref. [To] suggests that loops form at around 0.1 times the horizon scale, whereas Refs. [Tn[T2] 
(and some portion of the conventional wisdom) suggest a scale several orders of magnitude 
smaller. The larger loops would lead to enhanced signatures of several types, and tighter 
bounds on Gfi [T3] . 

In a system with a large hierarchy of scales, one might hope that analytic methods 
would be particularly useful in separating the processes occurring at different scales, while 
numerical methods would be needed to understand the nonlinearities at a given scale. This 
is the philosophy that we will attempt to implement. We focus on a microscopic description, 
similar in spirit to Ref. [H] and in particular Ref. [7] . An important difference from Ref. [7] is 
that we are less ambitious: that work was largely directed at understanding the horizon-scale 
structure quantitatively, whereas we are only interested in shorter scales. Also, we do not 
attempt to reduce the dynamics to a parameterized model but rather focus on the full two- 
point function; this two-point function appears to be characterized by a critical exponent, 
which did not enter into previous work. Finally, our (tentative) conclusions are opposite to 
Ref. [7], in that we believe that the network scales even without gravitational radiation. 

The full microscopic equations for the string network [TU [7] appear to be too complicated 
to solve, and so we model what we hope are the essential physical processes. We have had 

^For reviews of network properties and other aspects of cosmic strings see Refs. [31[S]. 
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some success in matching results from simulations, but there are also some discrepancies 
which may indicate additional processes that must be included. 

In section 2 we consider the evolution of a small segment on a long string. If the rate 
of string intercommutations is fixed in horizon units, then over a range of scales the only 
important effect would be stretching due to the expansion of the universe. We are then able 
to determine the two-point functions characterizing the small scale structure. We find that 
the string is actually rather smooth, in agreement with recent simulations [12]: its fractal 
dimension approaches one as we go to smaller scales. There is a nontrivial power law, but 
this appears in the approach of the fractal dimension to one; that is, the critical exponent x 
is related to the power spectum of perturbations on the long string. The two-point functions 
depend on two parameters that must be taken from simulations. One of these is the mean 
v"^ of an element of the long string; this determines the exponent x- The second parameter 
is the normalization of the two-point function. Our results for the two-point function match 
reasonably well with the simulations [12] over a range of scales, but there is substantial 
deviation at the shortest scales. It remains to be seen whether this is due to transient effects 
in the simulations, or real effects that we have omitted. 

In section 3 we study the effects of the small scale structure on string lensing. Because 
the string is fairly smooth, these effects are not as dramatic as has been considered in 
some previous work. However, there are calculable deviations from perfect lensing. We also 
consider, for lensing by a long string, the likely trajectory of the string relative to the axis 
of a given lensed object. We discuss some deviations from gaussianity in the small scale 
structure, due to kinks. 

In section 4 we add in loop production as a perturbation. Even though the strings are 
relatively smooth at short distance, we find that the total rate of loop production diverges 
for small loops; the rate of divergence appears to agree with recent simulations [11]. This 
divergence does not invalidate our stretching model, because the loop production is localized 
to regions where the left- and right-moving tangents are approximately equal, but it 
points to a complicated fragmentation process, cascading to smaller and smaller loops. We 
are not able to follow the fragmentation process analytically at present, but we give an 
analytic argument as to why the cascade should terminate, leaving loops at some small but 
fixed multiple of the horizon scale. We also point out an interesting puzzle related to loop 
velocities, which again points to a complicated fragmentation process. 

In section 5 we discuss various applications and future directions. The behavior of string 
networks is notoriously complicated, and it appears to remain so even when one focusses 
on small scales. Thus we view our work as part of an ongoing dialogue between analytic 
calculations and simulations, which we hope will lead to a more complete and precise picture. 
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2 The stretching regime 



2.1 Assumptions 

We consider "vanilla" cosmic strings, a single species of local string without superconducting 
or other extra internal degrees of freedom. The evolution of a network of such strings is 
dictated by three distinct processes. First, the expansion of the universe stretches the strings: 
on scales larger than the horizon scale irregularities on a string are just conformally 
amplified, but on smaller scales the string effectively straightens [15]. This is described by 
the Nambu action in curved spacetime. Second, gravitational radiation also has the effect of 
straightening the strings, but it is significant only below a length scale proportional to du 
and to a positive power of the dimensionless string coupling G/x [161 113 CHI [iHl Ej . Since 
this is parametrically small at small G/x, we will ignore this effect for the present purposes. 
Finally, intercommutations play an important role in reaching a scaling solution, in particular 
through the formation of closed loops of string. At first we shall neglect this effect but we 
will be forced to return to this issue in section 4. 

Let us consider the evolution of a small segment on a long string. We take the segment to 
be very short compared to the horizon scale, but long compared to the scale at which gravita- 
tional radiation is relevant. The scaling property of the network implies that the probability 
per Hubble time for this segment to be involved in a long string intercommutation event is 
proportional to its length divided by du-, and so for short segments the intercommutation 
rate per Hubble time will be small. Formation of a loop much larger than the segment might 
remove the entire segment from the long string, but this should have little correlation with 
the configuration of the segment itself, and so will not affect the probability distribution for 
the ensemble of short segments. Formation of loops at the size of the segment and smaller 
could affect this distribution, and the results of section 4 will indicate that the production of 
small loops is large, but this process takes place only in localized regions where the left- and 
right-moving tangents are approximately equal. Thus, there is a regime where stretching is 
the only relevant process. 

If we follow a segment forward in time, its length increases but certainly does so more 
slowly than the horizon scale dn-, which is proportional to the FRW time t. Thus the length 
divided by dn decreases, and therefore so does the rate of intercommutation. If we follow the 
segment backward in time, its length eventually begins to approach the horizon scale, and 
the probability becomes large that we encounter an intercommutation event. Our strategy 
is therefore clear. For the highly nonlinear processes near the horizon scale we must trust 
simulations. At a somewhat lower scale we can read off the various correlators describing the 
behavior of the string, and then evolve them forward in time using the Nambu action until 
we reach the gravitational radiation scale. The small probability of an intercommutation 
involving the short segment can be added as a perturbation. This approach is in the spirit 
of the renormalization group, though with long and short distances reversed. 
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2.2 Two-point functions at short distance 



In an FRW spacetime, 

ds'^ = -dt^ + a{tfdx ■ (ix 

= a{rf{-dr^ + rfx ■ dx) , (2.1^ 

the equation of motion governing the evolution of a cosmic string is [20] 



a e \ e 

Here e is given by 

X'2 XV2 



x+2-(l-x2)x=-f-V . (2.2) 



These equations hold in the transverse gauge, where x ■ x' = 0. Dots and primes refer to 
derivatives relative to the conformal time r and the spatial parameter a along the string, 
respectively. The evolution of the parameter e follows from equation (12. 2p , 

i = -2^i,\ (2.4) 
e a 

From the second derivative terms it follows that signals on the string propagate to the 
right and left with da = ±dT/e. Thus the structure on a short piece of string at a given 
time is a superposition of left- and right-moving segments, and it is these that we follow in 
time. In an expanding spacetime the left- and right-moving waves interact — they are not 
free as in flat spacetime. 

Prom Eq. (12.41) it follows that the time scale of variation of e is the Hubble time, and so 
to good approximation we can replace x^ with the time-averaged v"^ (bars will always refer 
to RMS averages), giving e oc a~^^ as a function of time onlyo From the definition of e it 
then follows that the energy of a segment of string of coordinate length a is E = fia{T)e{T)a. 
For simplicity we will refer to E/^ as the length / of a segment, 

/ = a{T)e{r)a , (2.5) 

though this is literally true only in the rest frame. The scale factor is 

a oc r oc r*^' , i^' = u/{l - u) , (2.6) 

where 

radiation domination: u = 1/2 , u = 1 , ^ 0.41 , 

matter domination: u = 2/3 , u' = 2 , v"^ 0.35 . (2-7) 



^Thc transverse gauge choice leaves a gauge freedom of time-independent a reparamcterizations. A 
convenient choice is to take e to be independent of a at the final time, and then e will be cr-indcpendcnt to 
good approximation on any horizon length scale in the past. 
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The mean RMS velocities for points on long strings are taken from simulations [12]. It 
follows that 

1(xt'^' (xt^ , C = (1 - "^v^y , C = (1 - 2t;^)z/ . (2.8) 

In the radiation era, ~ 0.1 and in the matter era Cm ~ 0.2. Thus the physical length of 
the segment grows in time, but more slowly than the comoving length [15], and much more 
slowly than the horizon length dn oc t. 

For illustration, consider a segment of length (10~^ to 10~'')(i//, as would be relevant for 
lensing at a separation of a few seconds and a redshift of the order of z 2^ 0.1. According to 
the discussion above, l/dn depends on time as t''"^ ~ t~^^^ in the matter era. Thus the length 
of the segment would have been around a hundredth of the horizon scale at the radiation-to- 
matter transition. In other words, it is the nonlinear horizon scale dynamics in the radiation 
epoch that produces the short-distance structure that is relevant for lensing today, in this 
model. This makes clear the limitation of simulations by themselves for studying the small 
scale structure on strings, as they are restricted to much smaller dynamical ranges. 

In terms of left- and right-moving unit vectors p± = x± ^x', the equation of motion (12.21) 
can be written as [1] 

P± T -P± = [Pt - (P+ ■ P-) P±] • (2.9) 

We will study the time evolution of the left-moving product p+(t, a) ■ p+(r, a'). For this it 
is useful to change variables from (r, a) to (r, s) where s is constant along the left-moving 
characteristics, s — s'/e = 0. Then 

^r(p+(s, r) ■ p+(s', r)) = -- ( p_(s, r) ■ p+(s', r) + p+(s, r) ■ p_(s', r) 

+a{s, r) p+(s, r) ■ p+(s', r) + a{s', r) p+(s, r) • p+(s', r) j , (2.10) 

where a = — p+ ■ p_ = 1 — 2t;^. 

The equations of motion (12.91 12.101) are nonlinear and do not admit an analytic solution, 
but they simplify when we focus on the small scale structure. If p_|.(s, r) were a smooth 
function on the unit sphere, we would have 1 — p+(s, r) • p+(s',r) = 0{[s — s']^) as s' 
approaches s. We are interested in any structure that is less smooth than this, meaning that 
it goes to zero more slowly than [s — s']^. For this purpose we can drop any term of order 
[s — s']^ or higher in the equation of motion (smooth terms of order s — s' cancel because the 
function is even). 

Consider the product p_(s,r) ■ p+(s',r). The right-moving characteristic through (s, r) 
and the left-moving characteristic through (s', r) meet at a point (s, t — 6) where 6 is of 
order s — s'Jfl Eq. (12. 9p states that p_|_ is slowly varying along left-moving characteristics 
(that is, the time scale of its variation is the FRW time t), and p_ is slowly varying along 

■^Explicitly, for given r we could choose coordinates where e(r) = 1 and s{T',a) — a + t' — t + 0{[t' — t]"^) , 
and then S = {s — s')/2. 
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right-moving characteristics. Thus we can approximate their product at nearby points by 
the local product where the two geodesies intersect, 



p_(s, r) ■ p+(s', r) = r - 5) + 0(3 - s') . (2.11) 

Then 

dr{p+{s,T) ■ p+(s',r)) = ^(a{s',r~ 6) + a{s,T + 6) 

-a{s, r) p+(s, r) ■ p+(s', r) - r) p+(s, r) ■ p+(s', r)) + 0([s - sf ) .(2.12) 

When we integrate over a scale of order the Hubble time, the 6 shifts in the arguments have 
a negligible effect 0{6) and so we ignore them. Defining 

h++{s, s', r) = 1 - p+(s, r) ■ p+(s', r) , (2.13) 

we have 



drh++{s,s',T) = --h++{s,s',T)[ais',T) + ais,T) + 0{[s-sf)] . (2.14) 
a 



Thus 



h++{s, s', n) = h++{s, s', To) exp|-z/' ' ^Hs', r) + «(s, r) + 0([s - sf )] | . (2.15) 

Averaging over an ensemble of segments, and integrating over many Hubble times (and 
therefore a rather large number of correlation times) the fluctuations in the exponent average 
out and we can replace a{s, r) with a = 1 — 2w^, 

(/i++(.,s',ri)) ^ (/i++(.,/,ro))(ri/ro)-2^''^ . (2.16) 

Note that in contrast to previous equations the approximation here is less controlled. We 
do not at present have a good means to estimate the error. It depends on the correlation 
between the small scale and large scale structure (the latter determines the distribution 
of a), and so would require an extension of our methods. We do expect that the error is 
numerically small; note that if we were to consider instead (In /i+_|_(s, s', r)) then the product 
in Eq. (12.151) would become a sum and the averaging would involve no approximation at all. 
Averaging over a translationally invariant ensemble of solutions, we have 

(/i++(a - a', n)) ^ (ri/ro)-^'^'^(/i++(a - a', ro)) . (2.17) 

We have used the fact that to good approximation (again in the sense of Eq. fl2.16p ). e is only 
a function of time, and so we can choose a = s — J dr/e and a — a' = s — s' . Equivalently, 

{hM^-cT',T))^l^^^. (2.18) 
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The ratio of the segment length to d/^ = (1 + u')t is 

/ ae(a — a') a — a' , ^ 

— oc ^ oc . 2.19 

The logic of our earlier discussion is that we use simulations to determine the value of 

at l/dn somewhat less than one, and then evolve to smaller scales using the Nambu action. 

That is, 

h++{a - a', r) = ho when a - a' = xot^^^'^'''' , (2.20) 

for some constants Xq and Hq. We assume scaling behavior near horizon scale, so that ho is 
independent of time. Using this as an initial condition for the solution fl2.18p gives 



(k,,(a-a',r)).h,[—^^ . Ml/ty^ ■ X=^^, = j—- (2.21) 

In the last form we have expressed the correlator in terms of physical quantities, the segment 
length / defined earlier and the FRW time 

Eq. fl2.2ip is our main result. Equivalently (and using a parity), 

(p+(a,r) ■ p+(a',r)) = (p_(^,r) • p^^',^)) ^ 1 - A{l/tf^ . (2.22) 

In the radiation era Xr ~ 0.10 and in the matter era Xm ~ 0.25. There can be no short 
distance structure in the correlator p+ ■ p_, because the left- and right-moving segments 
begin far separated, and the order a/ a interaction between them is too small to produce 
significant nonsmooth correlation. Thus, from (12. lip we get 

(p+(a, r) ■ p_(a', r)) = -a + 0{a - a') . (2.23) 



2.3 Small fluctuation approximation 

Before interpreting these results, let us present the derivation in a slightly different way. The 
exponent x is positive, so for points close together the vectors p+(cr, r) and p+((T',r) are 
nearly parallel. Thus we can write the structure on a small segment as a large term that is 
constant along the segment and a small fluctuation: 

p+(a, r) = P+(r) + w+(r, a) - ^P+(r)ti;2 (r, a) + . . . , (2.24) 

with P+(r)^ = 1 and P+(t) ■ w+(r, a) = 0. Inserting this into the equation of motion ( 12. 9p 
and expanding in powers of w+ gives 

P+ = -^[P_-(P+.P_)P+] , (2.25) 

W+ - = ^ [(w+ ■ P_) P+ + (P+ ■ P_) w+] 

= -(w+ ■ P+) P+ + -(P+ . P_) w+ . (2.26) 

(X 

"^Wc have not yet needed to specify numerical normalizations for cr and t, or equivalently for e and a. 
The value of ho depends on this choice, but the value of A does not. 
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Since the right-moving p_ is essentially constant during the period when it crosses the short 
left- moving segment, we have replaced it in the first line of ( 12.26^ with a a-independent 
P_. In the second line of fl2.26l) we have used w+ ■ P_|_ = 0. In the final equation for w_|_, 
the first term is simply a precession: P_|_ rotates around an axis perpendicular to both P+ 
and P_, and this term implies an equal rotation of w_|_ so as to keep perpendicular to 
P+. Eq. (12.261) then implies that in a coordinate system that rotates with P_|_, w+ is simply 
proportional to a~". 
It follows that 

(p+(a, r) ■ p+(a', r)) - 1 = -i ([w+(a, r) - w+(a', r)]') (2.27) 

scales as a~^" as found above (again we are approximating as in eq. (12.161) . and again this 
statement would be exact if we instead took the average of the logarithm). Similarly the four- 
point function of w+ scales as a"^"". We have not assumed that the field w+ is gaussian; the 
n-point functions, just like the two-point function, can be matched to simulations near the 
horizon scale. We can anticipate some degree of nongaussianity due to the kinked structure; 
we will discuss this further in section 3.2. 



2.4 Discussion 



Now let us discuss our results for the two-point functions. We can also write them as 

(x'(cr, r) ■ x'(cr, r)) 2[1 — v^j 

(x((T, r) • x(cr', r)) ^ /, , xo, / X 

corr,(/, t) ^ )} \ } '( ^ 1 - —,il/tf^ . (2.28) 
(x(cr, T) ■ x(cr, T)} Itr 

These are determined up to two parameters v"^ and A that must be obtained from simulations. 
A first observation is that these scale, they are functions only of the ratio of / to the horizon 
scale. This is simply a consequence of our assumptions that the horizon scale structure scales 
and that stretching is the only relevant effect at shorter scales. We emphasize that these 
results are for segments on long strings; we will discuss loops in Sec. 4. 

It is natural to characterize the distribution of long strings in terms of a fractal dimension. 
The RMS spatial distance between two points separated by coordinate distance a is 

(r2(/, t)) = (x' ■ x') f dl' f dl" corr,(/' - /", t) 
Jo Jo 

A{l/tf^ 



(1 - v^)f 



(2.29) 



(2x + l)(2x + 2)(l-t;2) 

We can then define the fractal dimension df (which is 1 for a straight line and 2 for a random 
walk) , 

2d In/ 



df 



d\n{r\l,t)) 
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Figure 1: Comparison of the model (dashed hne) with the data provided by [12] (sohd red 
hne) in the radiation-dominated era, for which the correlation length is ^ ~ 0.30t. 



The fractal dimension approaches 1 at small scales: the strings are rather smooth. There is 
a nontrivial scaling property, not in the fractal dimension but rather in the deviation of the 
string from straightness. 



We define the scaling dimension dg = 2%. Note that dg is not large, roughly 0.2 in the 
radiation era and 0.5 in the matter era, so the approach to smoothness is rather slow. 

Our general conclusions are in agreement with the recent simulations of Ref . [12] , in that 
the fractal dimension approaches 1 at short distance. To make a more detailed comparison 
it is useful to consider a log- log plot of 1 — corra;(/,t) versus /, as suggested by the scaling 
behavior f l2.3ip : we thank C. Martins for replotting the results of Ref. [12] in this form. The 
comparison is interesting. At scales larger than du (which is ~ 6.7.^ in the radiation era and 
~ 4.3^ in the matter era) the correlation goes to zero. Rather abruptly below du the slope 
changes and agrees reasonably well with our result. It is surprising to find agreement at 
such long scales where our approximations do not seem very precise. On the other hand, at 
shorter scales where our result should become more accurate, the model and the simulations 
diverge; this is especially clear at the shortest scales in the radiation-dominated era (Fig. [1]). 

One possible explanation for the discrepancy is transient behavior in the simulations. We 
have argued that the structure on the string is formed at the horizon scale and 'propagates' 




(2.30) 



1 -corr^(/,t) oc (//t)2^ , 1 



corrt(/,t) oc (l/tY'' . 
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Figure 2: Comparison of the model (dashed hne) with the data provided by [12] (sohd blue 
line) in the matter-dominated era, for which the correlation length is ^ ^ 0.69t. 



to smaller scales (in horizon units) as the universe expands. In Ref. [12] the horizon size 
increases by a factor of order 3, and so even if the horizon-scale structure forms essentially 
at once, the maximum length scale over which it can have propagated is 3^~^, less than half 
an order of magnitude. At smaller scales, the small scale structure seen numerically would 
be almost entirely determined by the initial conditions. On the other hand, the authors 
of Ref. [12] (private communication) argue that their result appears to be an attractor, 
independent of the initial conditions, and that loop production may be the dominant effect. 
Motivated by this we will examine loop production Sec. 4. Indeed, we will find that this 
is in some ways a large perturbation, but we are still unable to identify a mechanism that 
would produce the two-point function seen in the simulations. This is an important issue to 
be resolved in future work. 

Thus far we have discussed corr^,. Our result fl2.28p implies a linear relation between corr^. 
and corrt. In fact, this holds more generally from the argument that there is no short-distance 
correlation between p+ and p_, Eq. (12.231) : 



(1 — f^)corrj.(/, t) — v'^coTTt{l,t) = — 




P+{^, r) ■ P~{(t', r) + p_(a, r) • p+{a' , r) 



^ a . (2.32) 

Inspection of Fig. 2 of Ref. [12] indicates that this relation holds rather well at all scales 
below ^. 

The small scale structure on strings is sometimes parameterized in terms of an effective 
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tension j22| [23] . For a segment of length / the effective tension is given by 



where we have made use of result (I2.29P . Note that this is strongly dependent on the scale 
/ of the coarse-graining. 

In conclusion, let us emphasize the usefulness of the log-log plot of 1 — corr^,. In a plot 
of the fractal dimension, all the curves would approach one at short distance, though at 
slightly different rates. The difference is much more evident in Figs. 1 and 2, and gives a 
clear indication either of transient effects or of some physics omitted from the model. 



2.5 The matter-radiation transition 

In this subsection we will assume that our stretching model is actually valid down the scale 
where gravitational radiation sets in. If loop production or other relatively rapid processes 
are actually determining the small scale structure then this subsection is moot. 

We have noted that at very short scales we see structure that actually emerged from the 
horizon dynamics during the radiation era. Thus we should take the radiation-to-matter 
transition into account in our calculation of the small-scale structure. At the time of equal 
matter and radiation densities, 

1 - (p+(a, Teq) ■ p+(a', Teq)) ^ Mhjt.^f^^ , (2.34) 

where /gq is the length of the segment between a and a' at tgq. Assuming that the tran- 
sition from radiation-dominated to matter-dominated behavior is sharp (which is certainly 
an oversimplification), we evolve forward to today using the result (12.181) . The right-hand 
side of Eq. (12.341) is then multiplied by a factor (t/teq)"^*^™"™- In terms of the length today, 

I = 4q(V^eq)^") ^ave 

1 - (p+(a, r) ■ p^{a', t)) ^ Ml/tf ^^{t/t,^)-'^--''^^^-+'^^ . (2.35) 

This expression applies to scales l{t) that, evolved backward in time, reached the horizon 
scale dn before the transition occurred, i.e. at a time defined by /(t*) ~ (Ih such that 
< ^cq- For longer scales, which formed during the matter era (for which > tcq), we have 
simply 

1 - (p+(a, r) ■ p+(a', r)) ^ A^l/t)'''- • (2.36) 

The transition between the two forms occurs along the curve determined by the intersec- 
tion of the two surfaces (12.351) and (12.361) . This determines the critical length at the time of 
equal matter and radiation densities, /c(4q) = (^r/^m)^''^^'^'"~^'^"''*^eq- In terms of the length 
at some later time t, the transition occurs at 



0.002 



0123456789 10 

l/10"^d^ 

Figure 3: Structure on the string, as a function of the length / at present time (sohd 

curve). On scales larger than the critical length /<, ~ 3 x lO'^d^ the structure is determined 
by the matter era expression. On scales below Ic the transition result (12.351) gives an enhanced 
effect. The dashed curves show the extrapolations of the two relevant expressions: on small 
scales the actual structure is enhanced relative to the pure matter era result. 

so that the transition scale at the present time is /c ~ 3 x lO'^dn (Fig. [3]). The transition 
result fl2.35p implies more structure at the smallest scales than would be obtained from the 
matter era result, by a factor (/c//)^^^™"^'^ ~ (^c/0°'^- 

Of course, precise studies of the small scale structure must include also the effect of the 
recent transition to vacuum domination; this period has been brief so the effect should be 
rather small. 

3 Lensing 

Let us now consider the effect of the small scale structure on the images produced by a 
cosmic string lens. Previous work has discussed possible dramatic effects [2H , including 
multiple images and large distortions. We can anticipate that the rather smooth structure 
that we have found, which again we note is subject to our assumptions, will produce images 
with only small distortion. We will use our stretching model of the two-point function. If 
this proves incorrect one could apply the analysis using phenomenological values of x and 
A; for example, the extrapolation of the results of Ref. [I2] give a smoother string, and even 
less distorted images. 
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Figure 4: Schematic representation of the system considered, with the string lens displayed 
along the Xi-axis and the distant source and the observer located at points S and O respec- 
tively. 



3.1 Distortion of images 

We quote the result of ref. [26] for the angular deflection of a light ray by a string, 



7±(y±) = 4G'/i / da 



^''^7(0)7(0) x_L - 



X\\ 



|x± - y±| 



(3.i: 



J i=io{a) 



Here is the four-velocity of the unperturbed light ray, which we take to be (1, 0, 0, 1) as 
shown in Fig. HJ ^/'^ is a reference point on this ray, and subscripts _L and || are with respect 
to the spatial direction of the ray. Also, x'^{a,t) is the string coordinate|f| in terms of which 



Ft"" = x^x" - x^'x"' - -r]'"'{xPxp - xP'x'p) , 

2 



(3.2) 



and to{a) is defined by to(o") = ^^{a.tQ) — y^. 

In keeping with Sec. 2.3 we separate the string locally into a straight part and a fluc- 
tuation; we will keep the deflection only to first order in the fluctuation. We consider here 
only the simplest geometry, in which the straight string is perpendicular to the light ray 
and at rest, so that P+ = — P_ = (1,0,0). To first order in the fluctuation, x^{a,t) = 
{t, a, X2{cr,t) , X3{a,t)) . One then finds 



7i(y-L) = ^G/j, J da 
72(y±) = -sgn(?/2)/3 + 4^/2 / da 



-X3(o-, to){a - yi) + x'^ia, to)y2 



{a - yif + yl 

Xj,{a, to)y2 + X2{a, to){a - yi] 
(c^ - yif + yl 



(3.3) 



^ The region where the the light ray passes the string is small on a cosmological scale, so to use our earlier 
results we can locally set a — t~ 1, dl = da, dt — dr- 
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with [3 = AirGfi being half the deficit angle of the string. To the order that we work to is 
a constant, corresponding to the time when the light ray is perpendicular to the straight 
string. 

We can use our results for the small scale structure to calculate the two-point functions 
of the deflection. Let us focus on the local magnifications parallel and perpendicular to the 
string, given by basic lensing theory as 

Do oyi 

MM = 1 - ^'(^°-^') g2a(y,) . (3.4) 

Do oy2 

Here Do and Di are the distances of the source from the observer and the lens respectively 
(these would be the angular diameter distances on cosmological scales). It is particularly 
interesting to consider the differential magnifications for the two images produced by a string, 

SM, = Mi(y^) - Mi(yl) , 6M, = M^iy^) - M,{y'^) . (3.5) 

We take for simphcity y_L = (0,6) and y^ = (0, —b) for b = [3Di{Do — Di)/Do] this corre- 
sponds to the symmetric images of an object directly behind the string. Then 

From Sec. 2 we obtain 

(w+(o-,t) ■ w+(o-',t)) = (w_(o-,t) • w_(cr',t)) = 4 (x2(o-,t)x2(cr',t)) 

= A\{a/tf^ + {a'/tf^ - {[a - a']/tfA . (3.7) 



Note that Eq. (12.271) does not fully determine the two-point functions (13.71) . and in fact that 
the latter cannot be translation invariant. We have fixed the ambiguity by defining the 
expectation value to vanish when a or a' vanishes (that is, at the point on the string nearest 
to the light ray); this amounts to a choice of how one splits p± into P± and w±. 

As is stands, equation (13. 7p applies only when |cr — cr'| is larger than the critical length 
Ic- If this is not the case, one must take into account the radiation-to-matter transition as 
discussed in section 2.5. According to equation (I2.35p . the RHS of (13. 7p gets multiplied by 
a power of t/tcq- From {x^x^) we obtaiiu 

{6Ml) = {6Ml) = Ml:^MA^(26/to)2x.(t^/t^j-2C.n-2x.U+2x. (3,8) 

Plugging in the numeric values for x, A and C, and using a representative value for the 
dimensionless parameter ~ 10~^, we obtain a RMS differential magnification slightly 
below 1%: 

{5K'py'^ ~ 0.009 . (3.9) 



^The lensing scale, which is provided by 6, is typically of the order of 10~^(i//, thus much smaller than 
the critical length l^- We can then safely use the formula for {x'2x'2) valid for the smallest scales over the full 
range of integration. Corrections from the longest scales will increase the final result, but not significantly. 
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3.2 Alignment of lenses and nongaussianity 



Another question related to short- distance structure is the ahgnment of lenses. Suppose we 
see a lens due to a long string, with a certain alignment. Where should we look for additional 
lens candidates? Previous discussions [271 [28] have considered the two extreme cases of a 
string that is nearly straight, and a string that is a random walk on short scales; clearly the 
networks that we are considering are very close to the first case. 

We keep the frame of the previous section, with the lens at the origin in and aligned 
along the 1-axis. Then as we move along the string, the RMS transverse deviation is 



2 







= / da' da" {x2{a)x'2{a')) 



A]2 

{l/tf' . (3.10) 



4(X + 1 

The extension in the x direction is just /, so the RMS angular deviation is 

^A/A{x+l){l/t)^=I^ . (3.11) 

If we put in representative numbers, looking at an apparent separation on the scale of arc- 
minutes for a lens at a redshift of order 0.1, we obtain with the matter era parameters a 
deviation 5ip ~ 0.05 radians. That is, any additional lenses should be rather well aligned 
with the axis of the first. If the string is tilted by an angle ip to the line of sight, then 
projection effects increase 56 and by a factor 1/ cos^l). Of course, for lensing by a loop, 
the bending will be large at lengths comparable to the size of the loop. 

Lens alignment provides an interesting setting for discussing the nongaussianity of the 
structure on the string. If the fluctuations of x'2 were gaussian, then the probability of finding 
a second lens at an angle 5</9 to the axis of the first would be proportional to q-^'p'^ /"^^"p ^ and 
therefore very small at large angles. However, we have considered thus far a typical string 
segment, which undergoes only stretching. There will be a small fraction of segments that 
contain a large kink, and one might expect that it is these that dominate the tail of the 
distribution of bending angles. 

Let us work this out explicitly. Consider a left-moving segment of coordinate length cr, 
and let P(cr, r, k) dk be the probability that it contain a kink for which the discontinuity 
|p+ — p'_^| lies between k and k + dk < k < 2). There are two main contributions to 
the evolution of P. Intercommutations introduce kinks at a rate that we assume to scale, 
so that it is proportional to the world-sheet volume in horizon units, a^ear"'^'^' ^"dr, and 
to some unknown function g{k). Also, the expansion of the universe straightens the kink, 
A; oc a"" [IE]. Theifl 

1^.. ,3.12) 



^ Other effects arc often included in the discussion of kink density, such as the removal of kinky regions 
by loop formation but these have a small effect. 
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We set P to zero at the matching time tq defined by 

a = Xori+2-'^' , (3.13) 

as in Eq. (12.201) : earher kinks are treated as part of the typical distribution, while P identifies 
kinks that form later. For simplicity we assume that xq is small enough that the probability 
of more than one kink can be neglected. 
To solve this, define 

K = kr^' , Q((T, r, k) = kP{a, r, k) . (3.14) 

Then 

TdrQ = ag{KT~^')KT~^^''' . (3.15) 
This can now be integrated to give 



Jrn ^' 



a 



I TO 

dk' 



Here k^ = k{xo/aY'^^^~^^'^''"'^'' and a = a / t^^'^'^' '"^ . Note that a is just a constant times //t, 
so the probability distribution scales. The source g{k) vanishes by definition for /c > 2, so 
/co > 2 is equivalent to /cq = 2. 

Rather than the angle 5(p between the axis of the first lens and the position of the second 
lens, it is slightly simpler to consider the angle 69 between the two axes. In the small 
fiuctuation approximation this is just X2((j), and so the RMS fiuctuation is 

W = ^/A/2{l/tY . (3.17) 

In the gaussian approximation, the probability distribution is g-'^^^/^'^f and so is very small 
for large angles. On the other hand, a large angle might also arise from a segment that 
happens to contain a single recent kink. Treating the segment as straight on each side of the 
kink, the probability density is then precisely the function P(cr, r, k) dk just obtained, with 
k = 69. If we consider angles that are large compared to 69 but still small compared to 1, 
the range of integration in the solution ( 13.16^ extends essentially to the full range to 2 and 
so the integral gives a constant. Then 

P{a, r, k) dk oc k-^/^dk = 69-^/<d69 . (3.18) 

Thus the tail of the distribution is not gaussian but a power law, dominated by segments 
with a 'recent' kink. One finds the same for the distribution of 6(f. It is notable however that 
the exponent in the distribution is rather large, roughly 4 in the matter era and 10 in the 
radiation era. Thus the earlier conclusion that the string is rather straight still holds. The 
sharp falloff of the distribution also suggests that a gaussian model might work reasonably 
well, and indeed we will employ this in the next section. 
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4 Loop formation 



The model presented thus far assumes that stretching is the dominant mechanism governing 
the evolution of the small scale structure. Motivated by the discrepancy with Ref. [12] (and 
at the suggestion of those authors), we now consider the production of small loops. Given 
the smoothness of the strings on short distances as found above, one might have thought 
that loop production on small scales would be suppressed, but we will see that this is not 
the case. We will ignore loop reconnection, based on the standard argument that this is rare 
for small loops [2]. 

We will treat the stretching model as a leading approximation, and add in the loop 
production as a perturbation. However, we will find that the loop production diverges at 
small scales, and so the problem becomes nonlinear and our analytic approach breaks down. 
Ultimately we would expect that some sort of improved scaling argument can be used to 
follow the structure to still smaller scales, but this is beyond our current methods; it is 
necessary first to identify the relevant processes. 



4.1 The rate of loop production 

In this section we aim to determine the rate of loop production, which occurs when a string 
self-intersects. That is, we want to compute the number of loops dM with invariant length 
between / and / + dl originating from a self-intersection at a string coordinate between a and 
a + da during a time interval (t, t + dt). Loops much smaller than the horizon size evolve in 
an essentially flat space so we are able to use null coordinates u = t + a and v = t — a. In 
this case the left- and right-moving vectors p± depend only on n or f respectively and the 
condition for the formation of a loop of length I is L+ = L_ with 

ru+l 

L+(u,/) = / c/n'p+(u') , 

J u 

L_{vJ) = / dv'p.iv') . (4.1) 

Jv-l 

Hence the number of loops formed is given by 

dM = 6^^\L+{u,l) -L_{vJ)) \det J\dudvdl , (4.2) 

where J is the Jacobian for the transformation [u, v, I) L+ — L_. This formalism is as in 
Refs. [llllTj. 

Expanding around the stretching result in the form of section 2.3, we have the functional 
probability distribution 

^[P+,P-] ^ Po(P+,P-)n[w+]P-[w_] . (4.3) 

That is, with stretching alone the left- and right-movers cannot be correlated on small scales. 
We will further assume that V± are gaussian with variance given by the two-point function. 
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As discussed in Sec. 3.1 this two-point function is not completely determined by Eq. fl2.27p . 
and in general we have 



(w+(m') ■ w+K)) = -A\u' - u"\''^/e^ + f{u') + f{u") 
= 2G{u', u") . 



(4.4) 



The function G is defined for later reference, with the factor of 2 representing the sum over 
transverse directions. To specify the function / one needs to choose how to split p± into 
P± and w±. The natural choice here is to make P± parallel to L±, i.e. P± = L±/L±. The 
transverse fluctuations along the segment that forms the loop, say < a < I, then have 
vanishing average. 



u+l 



du' w I fn') 



dv'wJv') = . 



(4.5) 



Integrating Eq. (14.41) J^^^ du' and using the condition (14.51) to set this to zero gives f{u') 
(hiu' — u) where 



A 



cr 



2x+l 



+ (/-a 



(4.6) 



(2x + l)/t2x [" ■ - - 2x + 2 

Similarly, (w_(t;') • w^{v")) = 2G{v',v") with f{v') = (j){v - v'). 

The three-dimensional 5-function in Eq. (14.21) can be expressed in spherical coordinates 

L_) = -^^5(L+-L_)5(2)(P+-P_). (4.7) 
^+ 

We take the ^-direction to be aligned along P_|_ = P_, in which case the Jacobian matrix is 
given by 



as 



fwl{l)-wl{0) wl{l)-wl{0) l[wl{0)-wl{l)]\ 



\wl{l)-w^_{0) wl{l)-wl{0) l[wl{0) - wlil)] J 

For notational simplicity we have, after differentiating, set m = and v = I. The fluctuations 
w± effectively live in the x-y plane due to the orthogonality condition P± ■ w± = 0. 
Let us first estimate the scaling of dN' with /. We have 



L_ = - 
2 



da (wlia) - wl{a)) = 0{P^+^ /t'''') 



(4.9) 



The width of the distribution of — L_ is of order Z^"''^'*^ and so the average of (5(L+ — L_) 
is a density of order The (5-function (14. 7p then scales as while the Jacobian 

scales as l^^ (the columns are respectively of orders l^, l^, /^^) , giving in all 



dX 



oc 



dl 



IdM 



(X 



dl 



(4.10) 
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Thus, for 2x < 1 (as found in both the matter and radiation eras), the total rate of string 
loss diverges at small /: the calculation breaks down. 

We will discuss this cutoff further in the next subsection, but first we estimate the 
numerical coefficient. We see no reason for a strong correlation between the 5- function 
and Jacobian factors in dA/", so we take the product of their averages. For the radial part of 
the ^-function, a gaussian average gives 



dy^,y(^^L+-L-))^-4{{L+-L^)^)-iy4{(L+-L-f)+.. 
271 

^ P-J/2«(X)+0{J/^) 



e 

oo 27r 

^ (4.11) 



The subscripts c in the second line refer to the connected expectation values, obtained by 
contracting the gaussian fields w± with the propagator (14. 4p . We have defined 

R{X) = da da'G\a,a',t) = -^A'C{x) , (4-12) 

r(. ] - 1 f I 1 4 4r^(2+2x)\ 

^yX) - 4{l+2x)2 ll+4x^ {1+xr 3+4x r(4+4x) i " y^'^'^) 



Numerically C{x) = (0.0049,0.0106,0.0111) for x = (0.1,0.25,0.5). 

We now turn to the angular part of the probability distribution function. This must be 
a function of P+ ■ P-_, and we will take it to have exponential form [7]: 

P(P+, P„) = . ^ ^ e-^P+ P- . (4.14) 

(47r)2smhA 

The prefactor normalizes the distribution. The parameter A can be determined from the 
requirement that (P+ ■ P_) = —a , or equivalently 

1 cosh A _ /, 1 rN 

- - — — ~ -a . (4.15) 
A smh A 

Thus, A ~ 3a . More precisely, A is 0.55 in the radiation epoch and 0.95 in the matter epoch. 
We are now able to calculate the average value of (14. 7p . Recognizing that {L"^) — l"^ we 



have 



Ae"^ t^x 

(4.16) 



AtiA^AtxC{x) sinhA /^+^^ ' 
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We now consider (| det J|). In the gaussian approximation the relevant probabihty dis- 
tribution is 



P±(w±(/),w±(0)) = exp 



(2vr) 



where the index i is summed over the two coordinates x and y, and 

VIe 



<(0 

w^(0) 



(4.17) 



(4.18) 



(That is, tlie columns and rows of the 2x2 matrix M correspond to the points and /, not 
the index i.) As usual, the whole distribution is determined solely by the two-point functions 
which have already been determined in (14.41) and (14. 6p . M^^^*-' = {'V±iy±)'^)- Thus, one 
finds that 



Ail-x) \X 1 



det 



4x 



X 



1-x 



(4.19) 



We now have all we need to write down an expression for (| det J|) . Before we do so, let 
us perform a simplifying change of variables: 



Y± = 

under which the expectation value of the loop Jacobian takes a compact form: 



(4.20) 



det J|) 



2tt^ \t 



F(X) + i-^ G(X) 
1 + X 



(4.21) 



Here we have defined an 8-dimensional vector X = {X^, X_, F+, Y_, Z^, Z_, W^, and 
two functions of it: 



F(X) = (Y+W^ - Y^W+) {Y^ + Wl-Yl- Wl) , 
G(X) = 2{W+W. - F+F_)(X+ - X_)(Z+ - Z_) + 

+ {Y^W. + Y.W^) [(X+ - X.f - (Z+ - ZJf\ . 



(4.22) 



Numerically the integral is ~ 110 for the radiation epoch and ~ 94 for the matter epoch; 
analytic estimates agree. 

Finally, we can combine results (14. 16p and (I4.2ip . Noticing that dudv = 2dadt, we get 

[t 



d{M) 



da dt dl 



(4.23) 



with c = 0.121 in the radiation epoch and c = 0.042 in the matter epoch. 
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4.2 The small loop divergence and fragmentation 

The total rate of string loss J ld{Af) diverges at small / for x ^ 0.5, as is the case in both 
the radiation and matter eras. The fractal dimension approaches 1 at short distance, but 
the exponent x, which characterizes the approach to this limit, indicates a relatively large 
amount of short distance structure. For example, if we consider the functions p± on the 
unit sphere, then 6a oc 6p]!^: the effective fractal dimension is though the path is not 
continuous, there are gaps due to kinks. 

Of course, the total rate of string loss is bounded, and so for small loops we must be 
multiply counting. A precise account of this effect will be very complicated; see for exam- 
ple Ref. [7j, which parameterizes it in terms of a modification of the exponential distribu- 
tion fl4.14p taking the form of a 'hole' in the forward direction p+ = p_. For the present we 
just determine the scale at which our calculation must break down. In a scaling solution, 
the total amount of long string in a comoving volume scales clS ^oo ^ ^ 

/t"^, while stretching 

alone would give £oo oc a". The rate of string lost to loops must be 



hoops t 



Our result for the same quantity, with a cutoff / > /*, is 

dt) loops h, t^\t) {l-2x)t\t 

h ~ O.lSt (4.26) 



^ ' dl^l-] =. -^(-] . (4.25) 



Equating these gives 



in both eras. 

The cutoff (14.261) is rather large, just an order of magnitude below the horizon scale. 
However, this is just the beginning of the story. The formation of these large loops does not 
terminate the process of loop production on smaller scales; this continues on the large loop, 
unaffected at least initially by the separation of the loop from the infinite stringj§ Thus we 
must expect a rather dramatic fragmentation process. The divergence of the rate of loop 
loss implies that at scales / well below the naive number density of loops of size ~ / is 
greater than l/l: these become dense on the string. Thus we might come to the conclusion 
that the fragmentation process continues until we end up with all loops at the gravitational 
radiation scale, the opposite of the naive result fOej) above @ 



®By the same token, the fragmentation process begins even before the large loop pinches off. By causality 
the reconnection that forms the loop cannot immediately have an effect on the fragmentation process in its 
interior. Thus the large scale (|4.26p may not be seen in the primary loops as usually defined. Rather, we 
should follow null rays from the two world-sheet points of primary loop formation backwards to the point 
where they meet, and define the length of the 'causal primary' loop as twice the time difference so as to 
enclose the entire causally disconnected region. 

^The previous study [29] may provide little guidance. The small scale structure there was limited to 
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In fact, we believe that this is not precisely the situation either. Rather, a non-self- 
intersecting loop will occasionally form, likely an order of magnitude or more smaller than 
the primary loop, and this length of loop will then be lost to the fragmentation process on 
smaller scales. If a fraction l//c of the string is lost in this way at each scale, then the total 
length of string in loops of size less than / will scale as Z^/'^. The exponent may be rather 
small, so there could be a substantial production of small loops. The divergent rate of loop 
production implies that our initial attempt to use analytic methods to separate scales has 
been too simple, but we might hope that a more sophisticated scaling model along the lines 
just indicated will be successful. 

Let us briefly describe an analytic model of the formation of non-self-intersecting loops. 
When the loop forms, we can define the distributions 



Because of the high fractal dimensions of the curves p±(cr), we might think of these distri- 
butions as reasonably smooth functions. The structure on the high harmonics of the loop 
will tend to smear these distributions into gaussians, but they will be skewed by the partic- 
ular random values of the lowest harmonics. Thus the distributions p+ and p_ will differ. 
The production of the smallest loops is roughly local in p, and so it can only continue until 
whichever of the two densities is smaller has been depleted. This leaves residual distributions 



according to the sign of the difference at each point. This nonvanishing difference implies a 
cutoff on the production of the smallest loops, and allows non-self-intersecting loops of finite 
size to form. 

4.3 Comparison with simulations 

The results that we have found have some notable agreements, and disagreements, with 
simulations. One success is an apparent agreement with the recent simulations of Ref. [TT] 
for the distribution of loop sizes: that reference finds a number density of loops per volume 
and per length dn/dl oc with p = 2.5 ±0.1 in the matter era and p = 3.0 ±0.1 in the 
radiation era. We have exponents 3 — 2% = 2.5 and 2.8 respectively. Our exponent is for 
the production rate rather than the density, but in this regime the density is dominated by 
recently produced loops and so these are the same. 

Let us verify this, and also obtain the relative normalizations for the two quantities. The 
number of loops per comoving volume of length between / and I + dl is unaffected by the 

10 harmonics, and neither of the two spectra studied corresponds closely to our case: the less noisy was 
equivalent to x = 0.5, where the loop production is only logarithmically divergent, while the more noisy is 
dominated by harmonics near the cutoff and does not resemble our distribution. 




(4.27) 
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expansion of the universe, and we are assuming that we are at scales where gravitational 
radiation can be neglected. The number then changes only due to production: 

^ ^«^^l=^^f-V' , (4.29) 



dt\ di J 72^2 p yt 

where we have inserted the string scaling density 1/72^2. This integral is dominated by 
recent times as long as 3z/ — 1 — 2% > 0, as holds in both eras (again, a oc f^). Integrating 
over time then gives 

dn c / / \ ^^'''2^ 

(4.30) 



dl (3i^ - 1 - 2^)72^4 \t 
In the notation of Ref. ITTI. 



(3z/ - 1 - 2x)7 

where the last factor is from the conversion from dn to t. Using 7 = 0.59, 0.30 respectively 
in the matter and radiation eras [12], this gives Co = 1.25, 10.3. These are larger than found 
in Ref. [11] by factors of 20 and 260 respectively. A large difference is not surprising, as 
we have argued that there will be substantial fragmentation which will move the loops to 
smaller / where the distributions are larger. However, it remains to be seen whether this can 
account for such an enormous difference. The fragmentation need not change the scaling as 
long as the long-string correlation functions remains a power law, because the shapes of the 
produced loops, and the resulting fragmentation, all scale; thus the agreement of the loop 
production exponents may be real. 

As we have noted, the loop scaling must eventually break down due to conservation 
of string, but with these reduced normalizations this happens at much lower scales, ~ 
4 X 10~^t and 2 x 10~^t respectively. Below these scales we expect the distribution to 
fall as discussed in the fragmentation discussion. Thus we might expect the distribution 
to peak near these scales, but perhaps to be rather broad in both directions. A possible 
phenomenological formula for loop production after fragmentation would be 

d^^ c x2>^-3 ^ 

x = l/t . (4.32) 



dadtdl t3 1 + (x/x,)2x-i-i/fc 

The constant c is the reduced value of our c, determined by comparison with simulations such 
as ref. [11], the exponent 1/k, which is likely to be small, must be determined by improved 
simulations of fragmentation, and the value x* ~ is determined by conservation of total 
string. Specific signatures, however, might be dominated by the tails, either at large loops 
or small. 

The relatively large scale (14.261) is similar to that found in the recent paper [10]. As 
discussed in Sec. 4.2, we expect that the final loop size is much smaller because of extensive 
fragmentation; the divergence that we have found for small loop formation will be equally 
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present in the calculation of fragmentation of a large loop into smaller ones. Ref. [TO] finds a 
much less extensive fragmentation. This points up the need for more complete studies of this 
process. This should be possible within the context of flat spacetime simulations, beginning 
with long strings having the small scale structure as input (this does depend on resolving 
first the issue of the long string two-point function, to which we return below). 

The production of very small loops takes place when p+ ~ p_: that is, near cusps on the 
long strings as proposed in Refs. [30], [31] . Indeed, all the loops that we consider are rather 
smaller than the correlation length, so the functions p+ and p_ are each somewhat localized 
on the unit sphere, and necessarily in the same region since L_|_ = L_ . In this sense all of our 
loops are produced near cusps. It is possible that there is a second population of loops that 
form at scales much longer than the correlation length, though most of these will quickly 
reconnect with long strings. 

The picture of a complicated fragmentation process is consistent with the results of 
Ref. [12] . Now let us return to the two-point functions as found there. The partial agreement 
with our stretching model might now seem surprising, given the large production of small 
loops. If the loop formation were distributed on the long strings, there would be at least 
one additional effect, where the string shortens due to loop emission and more distant (and 
therefore less correlated) segments are brought closer together. However, the loop production 
is not distributed: it occurs when the functions p±(cr), as they wander on the unit sphere, 
come close together. Our picture is that whole segments of order (14.261) are then excised. 
Thus, most of the segments that remain on a long string at given time would have been little 
affected by loop production. 

The discrepancy at short scales, which is particularly evident in the radiation era (Fig. 1), 
remains a puzzle. Let us first note that at the shortest scales in Figs. 1 and 2 the slope of the 
curve approaches and possibly even exceeds unity, corresponding to the critical value x = 0.5: 
if the two-point function is of this form then the small-loop production converges. Thus it is 
appealing to assume that it is feedback from production of short loops that accounts for the 
break in the curve in Fig. iS This is distinctly different from the cutoff via a p+ — p_ hole 
as discussed in Sec. 4. 2 and and in Ref. [7]: the latter is a modification of the H — correlations 
(though one that would have little effect on the two-point function), whereas Fig. 1 shows a 
modification of the -|--|- correlation. 

We are unable to identify a physical mechanism that would account for the observed 
two-point function. Consider a very short segment on a long string. With time, its size 
as a fraction of the horizon size decreases, so it moves to the left on Figs. 1 and 2. If it 
experiences only stretching, it will follow the slope of the dashed line. The loop production 
must have a strong bias toward removing segments with large fluctuations to account for 
the simulations. However, the general picture above is that the initial loop formation occurs 
at the scale and so any small segments will be carried away without any bias on their 

""^^Note also that x = 0-5 corresponds to the functions p± being random walks on the unit sphere. This 
again suggests this value might be some dynamical fixed point. 
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own internal configuration. Thus we would expect the two-point function to be given by the 
stretching power law down to arbitrarily small scales, until gravitational radiation entersl"1 
We will not attempt to connect with the fiat spacetime simulations of Refs. [32l [10] 
because our approach requires expansion to drive the fractal dimension to 1 at short distance. 

Finally we would like to note a significant puzzle regarding the loop velocities. For a loop 
of length /, the mean velocity is 

V = ^(L+K /) + L.iv, I)) = yL+K /) . (4.33) 

Then 



III dada'{p^{a).p^{a')) 
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= 1 - 7 TT A^/t?'' ■ (4.34) 

(2x+l)(2x + 2)^/ ^ ^ ' 

Looking for example at loops with / = 10~^t, we obtain a typical velocity 0.985 in the matter 
era and 0.90 in the radiation era. These are significantly larger than are generally expected; 
Ref. [33] gives values around 0.75 and 0.81 for loops of this size. We emphasize that this is not 
a consequence of our dynamical assumptions, but can be seen directly by using the two-point 
functions of Ref. [12] (the discrepancy is then even greater in the radiation era). Assuming 
that both simulations are correct, we must conclude that the two-point function on loops is 
very different from that on long strings, in fact less correlated. One possible explanation is 
that loop formation is biased in this way, but this seems unlikely: loops form when right- 
and left-moving segments of equal L(/) meet, and this will happen most often in the center 
of the distribution. Rather, this seems to be another indication of the complicated nature of 
the fragmentation: the final non-self-intersecting loops must contain segments which began 
on the long string many times further separated than the final loop size. 



5 Conclusions and future directions 

Our attempt to isolate the physics on different scales has not been completely successful, but 
it suggests some interesting directions for future work. First of all, the two-point function on 
long strings must be better understood: is the discrepancy in Fig. 1 a result of omitted physics 
or numerical transients? This question can be settled numerically, by studying sensitivity to 
initial conditions, and by studying the time-dependence of an ensemble of small segments. 
A good understanding of the two-point function would then allow extrapolation to much 

^^Onc effect that may increase x slightly is that segments that do not fall into the loop production 'hole' at 
p+ ~ p_ may experience on average a larger value of a = — p+ -p- and so have a reduced two-point function 
— see Eq. (|2.15p . However, it seems impossible that this could account for the factor of 5 discrepancy of 
slopes evident in Fig. 1. 
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smaller scales, for application to lensing (as in section 3), gravitational radiation [10], and 
possible interference between short distance structure and cusps [3T] . 

Second, our work points to correlations among the two-point function, loop production, 
fragmentation, and loop velocity. In particular, our conjecture of a very complicated frag- 
mentation process can be tested — for example, the argument that non-self-intersecting loops 
must contain segments that originated far apart on the long string. Also, the normalization 
difference between our loop production and that in Ref. [H] must presumably be understood 
in terms of fragmentation. 

In conclusion, cosmic string networks are probably not as complicated as turbulence, but 
they share the property that a rather simple set of classical equations leads to a complicated 
dynamics that challenges both numerical and analytical attacks. We hope that our work is 
a step toward a more unified understanding of the properties of these networks. 
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